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I. INTRODUCTION 

A possible way to prepare ultra-cold atoms in strongly 
correlated states is to use magnetic Feshbach resonances 
P, Q • Using this techniques in Bose gases a breakthrough 
has been achieved with the observation of three-body Efi- 
mov states and also four-body resonances tightly 

linked to the Efimov effect Q- However, a stable 
many-body state has never been observed for bosons at 
resonance as a consequence of a large losses into deeply 
bound states [sj. In all these experiments, atoms are neu- 
tral and the pairwise interaction between them is char- 
acterized by a van der Waals tail. In the vicinity of a 
magnetic Feshbach resonance, the scattering length (de- 
noted by a) can be tuned (formally from —oo to oo) via 
an external magnetic field of amplitude denoted by S. In 
the vicinity of the resonance it can be parameterized as 



(1) 



where Bq is the position of the resonance, /S.B is the 
width of the resonance such that a = for i? = ;Bo + AS 
and flbg is the background scattering length, i.e. the 
off-resonance s-wave scattering length. Away from reso- 
nances the two-body scattering length is of the order of 
the van der Waals range 
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where m is the mass of one atom and Cq is the London's 
constant. In the resonant regime, this length is small as 
compared to the absolute value of the scattering length 
-Rvdw ^ |a|. The van der Waals length gives the short 
distance scale for interatomic forces and also permits to 
define a low energy regime corresponding to energies E 
such that |£^| < (h^/mRl^^). 



The precise description of the interaction between ul- 
tracold atoms needs multi-channel computations. In 
this paper we use instead a separable two-channel model 
which describes the Feshbach mechanism and contains 
the relevant energy scales in the vicinity of the resonance 
while short range properties are taken into account only 
qualitatively. It is a simplified version of the model of 
Refs. [H, m and it was already introduced in the con- 
text of the two spin-component interacting Fermi gas and 
for the study of the three- and four-body problems in 
Refs. hj-nd. A similar version of the model has been 
also used in the context of Efimov physics for three spin- 
component fermions |17| . It has also the same structure 
as the one introduced in the study of three interacting po- 
larized fermions in Ref. [18]. In Refs. it has been 
shown that this separable two-channel model allows one 
to study the few-body problem by solving equations of 
the same order of difficulty than in the zero range ap- 
proximation. The model takes into account the direct 
interaction between atoms and is thus able to describe 
the background scattering and the finite width of the 
resonance. In this paper, we use this feature in order to 
analyze the respective contributions of the background 
scattering length, of the width of the resonance and of 
the finite potential range on the properties of the shal- 
low dimers in the different regimes encountered in experi- 
ments using the magnetic Feshbach resonance techniques: 
broad and narrow resonances and also the possible neigh- 
borhood of shape resonance. Deviations from universal 
laws valid near the threshold of resonance are obtained in 
these three regimes. We use the example of a broad res- 
onance in potassium where the spectrum has been com- 
puted in a microscopic coUisional approach in Ref. |19j to 
show that the basic description of the short range interac- 
tion with a single parameter is enough for a quantitative 
description of low energy properties. The case of narrow 
resonances, is particularly interesting because these de- 
viations can be evaluated accurately without explicitly 
taking into account short range physics. For this type of 
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resonances, we define an intermediate regime where the 
binding energy of the dimer is an affine law of the exter- 
nal magnetic field. Finally, we show that the two-channel 
model permits to study narrow resonances for a larger in- 
terval of magnetic detuning as compared to the effective 
range approach pO | . Finally, we find an approximate law 
for the recombination constant which takes into account 
the background scattering length Obg and also general- 
izes the analogous law obtained in the framework of the 
effective range approximation. 

The paper is decomposed as follows. In section HIl we 
define the two-channel separable model used in this pa- 
per. In section Hill we analyze the basic properties of the 
shallow dimers in the different regime of the the magnetic 
Feshbach resonance. In section ITVl we derive the integral 
equation which permits to solve the three-body problem. 
We then consider the recombination rate issue for three 
identical bosons and apply the formalism to the case of 
the intermediate regime of narrow resonances defined in 
this paper. 



II. TWO-CHANNEL SEPARABLE MODEL 

We consider neutral and spinless bosonic atoms. They 
belong to the open channel of the model characterized 
by a continuum of delocalized states above zero energy. 
The model supports a closed channel characterized by 
a discrete and structureless bosonic molecular state. In 
our terminology, molecules in the closed channel are dis- 
tinct from the dimers which are the bound states of the 
full Hamiltonian. The Feshbach resonance mechanism 
is encapsulated by a coherent coupling between atomic 
pairs of the open channel and the molecular state. The 
model contains a direct interaction term between atoms 
in the open channel which permits to take into account 
the background i.e. off- resonance scattering length (de- 
noted by ftbg). 

We use the second quantized form, where the oper- 
ator Ok annihilates one atomic boson of mass m and 
wavevector k in the open channel while 6k annihilates 
one molecule of mass 2m and wavevector k in the closed 
channel. Both operators at and 6k obey standard bosonic 
commutation rules: 



[a^^aU = (27r)'^J(k-k') , [6k, 6^1 = (2^)'<5(k-k'), (3) 

where the factor (27r)'^ in front of the (5-distribution cor- 
responds to the choice (r|k) = exp(ik • r) for the plane 
wave. Any other commutator vanishes. In what follows, 
the kinetic energy of an atomic plane wave of wave vector 
k is 
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The Hamiltonian of the model can be written as the sum 
of three terms: 



H = Hat + Hjjioi + Hi 



at — mob 
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where i?at is the Hamiltonian for bosons in the open 
channel, i?nioi is the Hamiltonian for molecules and 
-ffat-moi couples any atomic pair with the molecular 
state. The three terms in Eq. ([S]) are defined as follows: 

• The Hamiltonian iJat describes atoms: 
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For simplicity, the pairwise interaction among 
atoms has been chosen separable: it is character- 
ized by a coupling constant and a normalized 
function Xk with a Gaussian shape and a cut-off 
length denoted by 6: 



Xk = exp(-fc^6V2) 



(7) 



The parameter 6 mimics the potential range of the 
true interaction between neutral atoms and is thus 
of the order of the van der Waals range R^dw- The 
particular choice for the shape of the function Xk is 
not essential: it allows to perform some analytical 
simplifications and gives a qualitative description 
of the short range physics. The model Hamiltonian 
Hat allows one to describe two-body scattering far 
from the Feshbach resonance: the background scat- 
tering length of the model is a function of go and 
6. 

• Molecules in the closed channel are described by 
-ffmoi: 
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where i?moi is the internal energy of a molecule de- 
fined with respect to the zero energy of the open 
channel. In the magnetic Feshbach resonances, 
En\o\ is tuned by use of the external magnetic field 
B. We make the reasonable assumption that in the 
vicinity of the resonance, it can be approximated 
by an affine law: 



(9) 



In Eq.®, Bf is the magnetic field at which the 
molecular state energy crosses the threshold of the 
open channel continuum and denotes the differ- 
ence between the magnetic moments for an atomic 
pair in the open channel and the molecular state in 
the closed channel. 

• The last term in Eq. ([S]) couples the two channels 
and thus models the Feshbach resonance mecha- 
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In Eq. pOj) . the coupling function is generic: it 
is characterized by a constant A (chosen here to 
be real) and the same cut-ofF function Xk as in 
Eq. (O. When applied in the context of the few- 
body problem, this particular choice greatly sim- 
plifies the original model of Ref. it permits 
to obtain single closed equations for the three- and 
four-body problems [16, 18] similar to the ones ob- 
tained in the zero-range potential approach [21j]. 



III. TWO-BODY PROPERTIES 

A. Scattering amplitude 

The basic object for the study of two-body proper- 
ties is the scattering amplitude. Thus, in this subsection 
we derive the two-body scattering amplitude of the two- 
channel model to make the link between two-body prop- 
erties which may be measured experimentally {a,a\^g,AB, 
Edim) and the four parameters of the model (170, A, 

The most general state for the two-body system in the 
center of mass frame is a coherent superposition of two 
atoms in the open channel and of one molecule at rest in 
the closed channel: 
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^(k) is thus the wavefunction for a pair of atoms. The 
equations verified by the amplitudes ^(k) and /3 in 
Eq. (jlip are obtained from the projection onto the atomic 
and molecular subspaces of the stationary Schrodinger's 
equation at energy E: {E — H)\^) — 0: 

{E - 2ek) Ak = Axk/3 + 5oXk7 (12) 
{E - En,,i)l3 = 2A7, (13) 

where 7 = / ^^Xk'^k'- From Eqs. (|12ll3p we elimi- 
nate 7 and obtain: 

[E - 2ek) Ak~^(E- E^oi + —] /? = 0. (14) 
2A V 90 J 

For a scattering process, the atomic wavefunction is char- 
acterized by an incoming plane wave of momentum ko 
and an outgoing spherical wave: 



Ak = (27r)3<5(k - ko 
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where we use the standard prescription E E + iO^ to 
ensure that the scattered wave is indeed outgoing. Gen- 
eral scattering theory relates the scattering state j^") 
with the incoming state l^-^^)) by \^) = {1 + GoT)\^^°^), 
where T denotes the transition operator and Gq is the 
resolvent of the non- interacting Hamiltonian [23]. From 



this identity and from Eq. (|15l) . the half on-shell transi- 
tion matrix is: 

{k\T{E + zO+)]ko) = ^(e- E„,oi + —) p. (16) 
2A V 50 / 

The scattering amplitude f{E) is derived from the on- 
shell transition matrix (fc = fco) using the relation: 
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This way, f{E) can be expressed as 
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Equation (IT8|) shows that in this model scattering only 
occurs in the s-wave channel which is a sufficient descrip- 
tion of binary collisions in the regime of temperature and 
densities where ultracold gases are studied. In the ab- 
sence of the inter-channel coupling (A = 0) i.e. in the off- 
resonant regime, the scattering amplitude at zero energy 
is by definition equals to the opposite of the background 
scattering length Obg and from Eq. (|18p . one finds: 
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The system is in the neighborhood of a shape resonance 
when ]abg] ^ b and in this model this corresponds to the 
situation where go ^ .g§. Away from a shape resonance, 
]abg] is of the order of the potential range b. 

The general calculation of the scattering amplitude at 
finite energy in Eq. (|18p is easy to perform in the domain 
of negative energy E — h^kQ/m < 0. Keeping in mind 
the standard analytical continuation ko — iq with q > 0, 
one obtains: 
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(20) 

where erfc() is the complementary error function which 
equals unity at zero. The zero energy limit of Eq. (j20p is 
equal to the opposite of the inverse scattering length that 
can be exactly identified with the standard phenomeno- 
logical expression of Eq. ([T]). This provides the expres- 
sion of the width of the resonance as a function of the 
parameter of the model: 
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and also the magnetic detuning from the center of the 
resonance: 
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In what follows, Sfj.{B — Bq) is called the energy detuning. 
Equation (l?T|) shows that SfiAB has the same sign as the 
background scattering length Obg- We have checked that 
this property is indeed verified for the various resonances 
reported in Ref . [ISJ . As a consequence of the interaction 
between the open and the closed channel the resonance 
is not located at Sq' and from Eqs. (|9l22p the resonance 
shift in the magnetic field is: 



AB 
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One can notice that nearby a shape resonance where 
|oibg| ^ b the resonance shift appears as independent 
of the short range details of the interatomic forces 
and equals —AB. This is the case of the the Fesh- 
bach resonance in cesium at —12 G where one finds 
{Bq — Bo)/ AB ~ —1.12 if one takes the reasonable value 
b = i?vdw for the potential range. 



B. Dimers 

1. Bmdmg energies 

Using the results of the previous section, we now dis- 
cuss the properties of dimers i.e. the two-body bound 
states of the full Hamiltonian in Eq. ([S]). Their binding 



energies (denoted by E'dii 



> 0) are given by the 



poles of the scattering amplitude f{E) at negative energy 

{E = -Edira < 0) I 

following equation: 



{E = -Sdim < 0) [221. From Eq. ^ this leads to the 



gdimflbg e'^dim'' erfc(gdim&) 
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In the limit where E and S^{B — Bq) have both large and 
negative values, the leading term in the expansion of the 
left hand side gives the solution E — Enwi = ^-^dim: *-e., 
far from resonance the dimer coincides with the molecular 
state and is essentially in the closed channel. At finite 
detunings, one may distinguish two qualitatively different 
regimes in Eq. ([M|) : 

(i) For Obg < 6\/7r, the left hand side of Eq. (|24| is 
negative and is a monotonically decreasing func- 
tion of qdim- Hence, there is at most one possi- 
ble solution which only exists for a negative energy 
detuning d^j,{B — Bq). The dimer results from the 
inter-channel coupling and is denoted below as the 
Feshbach dimer. 

(ii) For Obg > 60r, the left hand side of Eq. ([24]) di- 
verges at a given value denoted in the following 
by g'*''^. It decreases from to — oo in the in- 
terval Qdim € [0,q'*'^[ and from +oo to —oo in the 



interval gdim €]g ^^,-I-cxd[. Consequently, there al- 
ways exists a solution in the interval -|-oo[. 
Away from the Feshbach resonance and for posi- 
tive energy detunings [S^{B — Bq) — >■ -l-oo], this so- 
lution gives the energy of the dimer associated 
with the direct coupling and we denote it as the 
background dimer. For negative energy detunings 
[5fi{B ~ Bq) < 0] a, second dimer exists: at thresh- 
old [i.e. for SiJ,{B — Bq) — 0^] it corresponds to the 
Feshbach dimer and for decreasing values of the en- 
ergy detuning there is an avoided crossing between 
the Feshbach and the background dimers. 

Due to their important role and their universal character, 
we now focus on the low energy two-body properties of 
this model. We thus assume that the dimer's energy is 
such that Qdimb ^ 1. In order to obtain an approximation 
beyond the universal regime where gdim — 1/a, we solve 
the equation: 



E + S^ijBQ -B + AB) 
W) 



0, 



(25) 



by expanding expanding the left hand side up to second 
order in q and find: 
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(26) 



where the parameter R* denoted hereafter as the width 
radius is linked to the width of the resonance by : 



R* = 



mahgSfj.AB 



(27) 



We anticipate that the expression in Eq. (|26|) cannot be 
obtained from the effective range approach (even in the 
formal limit where 6-^0). In the limit of large and pos- 
itive values of the scattering length, i.e. at the thresh- 
old of appearance of the Feshbach dimer, Eq. (pS)) gives 
at the lowest order the universal and standard law sup- 
ported also by the Bethe-Peierls zero-range model [l^j : 
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(28) 



The interest of the two-channel model is to give quan- 
titative corrections to this law. In the vicinity of the 
threshold where a is much larger than the three length 6, 
R* and |abg|, the correction in 1/a^ on qdim is given by: 



1 



gdii 



1 



a 



2b 



(29) 



where terms of higher order are not consistent with the 
approximation made for Eq. (|26p . A similar law was de- 
rived by using another two-channel model in Ref. (isj 
(see Eq. (53) in this last reference). We now distinguish 
different regimes where the deviation to the universal pre- 
diction are qualitatively different. 
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a. Broad resonances- For a broad resonance, R* is 
of the same order of magnitude or smaher than 6, mean- 
ing that the correction to the universal law depends on 
the details of the true interatomic potential and Eq. 
is valid in the regime of small detuning: 



\B-Bo\ 
\AB\ 



< 1. 



(30) 



far from the Feshbach threshold the background dimer is 
shallow. In the regime where a of the same order of 
magnitude than Obg, the expansion of Eq. (|26p for large 
values of abg (for a fixed ratio Ohg/a) gives the correction 
to the universal law gdim ^ l/obg: 
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b. Narrow resonances- The case of narrow reso- 
nances defined by R* ^ |abg| has attracted some interest 
in the three-bosons problem [2^ [2^ and we focus on 
the dimer spectrum given by the present realistic model. 
Near threshold for this type of resonance, the regime of 
validity of Eq. ([25)1 is very tiny: 



\B-Bo 
\AB\ 



|Qbg| 
R* 



< 1. 



(31) 



However for a narrow resonance, the first correction to 
Eq. (pS)) given by Eq. (PU)) has also an universal character: 
a feature already noticed in Refs. [H, H^]. A more inter- 
esting feature brought by the two-channel model concerns 
higher relative detuning \B — So|/|AS|, where there still 
exists a shallow dimer. Equation (|26p permits to distin- 
guish between the contributions of the background scat- 
tering length flbg and of the potential range b. In the 
regime where R* 3> b | | , terms involving b are neg- 
ligible and one can thus use the following approximate 
law: 



,+ + 4i?*(a- abg) 



the 



Moreover, 
R*{a — flbg) » the 
bach dimer is given by: 



2i?*(a - flbg) 
intermediate 



regime 
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dim 



Sfi{B„ - B). (33) 



This law, derived in the limit of a large width radius R* , 
shows that for a negative background scattering length 
(obg < 0) and for a vanishing value of the scattering 
length (|a| <C |abg|), there still exists a shallow dimer of 
binding energy i?dim — —SfiAB. In terms of magnetic 
field detuning, the definition of the intermediate regime 
chosen here can be written as: 



R* 

I abg 
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It corresponds to a larger interval of detuning than the 
intermediate regime of Ref. [l^l where it is defined as 
^ jobgl and corresponds to the inequalities: 



R" > \a 

> 



AB 



B-Bo 



> 1. 



c. Vicinity of a shape resonance- In the case where 
the system is in the vicinity of a shape resonance i.e. 
|ctbg| ^ b and for a positive background scattering length, 



The presence of the parameter b in Eq. ([55)1 shows that 
this correction depends on the details of the real inter- 
atomic forces at short distances. For arbitrary magnetic 
detuning, the set of solutions of Eq. p4| can be ob- 
tained by using a graphical reasoning. As an example, in 
Fig. ([IJ we have plotted the possible values of gdimObg in 
Eq. (p4)) as a function of the energy detuning S^{B — Bq) 
for a specific value of the ratio \at,g\/R*. For a negative 
background scattering length (abg ^ ~b) only Feshbach 
dimers exist and the corresponding branch is plotted in 
the negative region of the product gdimflbg and of the 
energy detuning. For a positive background scattering 
length (abg ^ b), gdimflbg is positive and there are two 
possible branches displayed in Fig. ([T]). In this regime, 
the figure nicely illustrates the interplay between the Fes- 
hbach dimer and the background dimer [2^ . At the Fes- 
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FIG. 1: Graphical solution of Eq. (|24|) for Qdimb ^ 1 and the 
ratio |abg|/-R* = 2. The product gdimibg is plotted as a func- 
tion of the energy detuning (5/i(B — So). In the regime of neg- 
ative background scattering length (abg < 0) only one shal- 
low dimer exists (solid line). For positive scattering length 
(flbg < 0) there exists an avoided crossing between the Fesh- 
bach dimer and the background dimer respectively (dashed 
and dotted lines respectively). 

hbach threshold \&[i{B — So) = 0], the binding energy of 
the dimer of finite energy is: 



-Edim 



4TOaL 



1 + Jl + 4abg/i?* 



(36) 



Hence, as we consider the vicinity of a shape resonance 
where |abg| ^ Rvdw, this bound state is shallow only for 
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a sufBciently narrow resonance. This shows the possible 
occurrence of a regime where two low energy dimers co- 
exist. Interestingly, in the case of a narrow resonance in 
the vicinity of a shape resonance Eq. (1321) is a low energy 
solution for all values of the magnetic detuning including 
the off- resonant regime where a is of the order of abg- 



2. Occupation probability in the closed channel 

The inter-channel coupling induces a non vanishing oc- 
cupation probability of the dimer in the closed channel. 
In what follows, this probability is denoted by Pciosod- 
For a normalized bound state: 



(37) 



and Pciosed — In the case of a shallow dimer such 

that Qdimb <C 1, this probability is independent of the 
short-range parameter b: 



Pcloscd 
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(38) 



Near the Feshbach threshold, for a vanishing and nega- 
tive value of the energy detuning 6fi{B — Bq), the Fesh- 
bach dimer has a vanishing occupation probability in the 
closed channel: 



2R* 



Pcloscd ■ 



(39) 



while in the intermediary regime of Eq. p4p . this occu- 
pation probability is close to unity: 



a broad resonance, where the dimer near threshold is 
essentially in the open channel [l^ [l^l . Comparing this 
law to Eq. (j29p confirms our assumption on the parameter 
b. In the case of narrow resonances where R* is large as 
compared to |abg|, Eq. p6| shows that low energy two- 
body properties depend only sightly on the precise choice 
made for b which is of the order of R^dw- 

For a precise determination of 6, one may use data 
provided by detailed collisional models or experimental 
results when available. As an example, the spectrum for 
the resonance at 402 G in '^^K has been computed in 
Ref. [l^ (see Fig. 4 in this last reference). It appears 
that the choice b = 1.2i?vdw permits to fit precisely the 
spectrum even for a relatively large detuning. We plot- 
ted in Fig. ([2]) the different predictions: the semi-classical 
law (HTI) . the universal law (pS)) . the microscopic compu- 
tation of Rcf. [19] and the result given by the present 
two-channel model. In Ref. UM , we have also considered 




1 / Q - Qbg 
Pclosod ^ ~ 2 V R* ' 



(40) 



C. Short range parameter of the model 

In standard broad resonances used in current experi- 
ments, for a detuning |B — So| of few Gauss and even if 
-E'dim < ^-vdWi the deviation of the dimer's binding en- 
ergy from the universal threshold law of Eq. (j28p is non 
negligible. Interestingly, this property can be used as a 
way to set a reasonable value for the short range param- 
eter b. In order to justify that b is of the same order of 
magnitude than i?vdWj we use an approximate formula 
which has been derived semi-classically for a one-channel 
model with a potential containing a van der Waals tail 



FIG. 2: Spectrum of the Feshbach dimer for the resonance 
at 402 G in ^^K. Stars: microscopic computations [lH; con- 
tinuous line: two-channel model used in this paper where 
b = 1.2 f?vdw; dashed line: universal law i?dim = /{ma^)\ 
dashed dotted hne: semi-classical prediction from Eq. (|41|) . 
Inset: for large detunings, the binding energy tends to — -Bmoi- 



the resonance at -12 G for ^'^'^Cs where experimental re- 
sults of spectroscopy were obtained [131 • We have shown 
also that the data can be reproduced by a specific value 
of the short range parameter which is close to the van 
der Waals range (6 = 0.7 -Rvdw)- These studies provides 
strong confidences that for the various magnetic Fesh- 
bach resonances that can be encountered in experiments, 
the model can reproduce accurately two-body properties 
for energies lesser than E^dw- 



_ r(3/4)i?vdw /.-.N 
where, a= — ^ . 41 

V2r(5/4) 



It has been shown that this approximation is relevant 
in the limit of large scattering length a 3> Rvdw a-nd for 



D. Effective range approximation 

In the effective range approach, the inverse scattering 
amplitude 1/ f{E) is replaced by its approximation at the 
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linear order in E: 



for E <0. 



(42) 



In Eq. To is the effective range and q = —ik ii E > 0. 
In what follows, we compare the predictions for the shal- 
low dimer obtained from Eq. (|42|) to the one given by 
the two-channel model, we thus choose the value of the 
effective range obtained from Eq. (|20l) : 
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7^ 



262 
a 



■0 



2R* 1 



abgV 
a / 



(43) 



At resonance the effective range is r™*' = — 2R*, thus 
Tg®'' has a positive upper-bound of the order of the short 
range length 6 but depending on the parameters of the 
model, it can in principle assume any negative value. For 
a broad Feshbach resonance the effective range is of the 
order of the van der Waals range while for a narrow Fesh- 
bach resonance {i.e. R* \ahg\) it is large and negative 
^ —2R*) with a value which only depends on the 
general characteristics of the resonance through Eq. p7| . 

In the regime where the scattering length is large and 
positive, the effective range approach supports a shallow 
bound states with a binding wavenumber given by 



9dim — 



- Va^ - 2c 



(44) 



Equation (HH) and Eq. (1^51) coincide the each to the 
other only in the region of asymptotically small detunings 
\{B — Bq) / IS.B\ <C 1 of narrow resonances (characterized 
by a large width radius R* ^ b). In this very specific 
regime, the binding wavenumber of the Feshbach shallow 
dimer is 



9dii 



. + Va^+ 4:R*c 
2R*a 



(45) 



which differs from the more general approximate law in 
Eq. ([5^ . The effective range approach is thus especially 
relevant for narrow resonances and for |a| 3> |abg|, while 
for broad resonances it can only supports a qualitative 
description of the corrections to the universal laws in 
Eq. 

One can also compare Eq. ((24| for the binding 
wavenumber, with the equation written in a similar form 
but obtained in the framework of the effective range ap- 
proximation. In the low energy limit where qdimb <Si 1 
and in the effective range approximation: 



(Zdimflbg 



"■ ydim 



1 



<?dimabg - 1 md^AB (l 



B-B, 



AS") (1 " 9dimabg) 

^ B~Bo 
AB ■ 



(46) 



Equation (j46p coincides with Eq. (1^^ only in 
the regime of asymptotically small detuning i.e. 



|(i3-Bo)/AS| < 1) and for gdim|abg| < 1- Conse- 
quently, the interplay between a shape and a Feshbach 
resonance cannot be properly described by a model rely- 
ing on an effective range approach. 

In Fig. ^ we have plotted the dimer spectrum for the 
narrow resonance of "^^K at 752 G obtained from the two- 
channel model. We have added the different approximate 
laws of Eqs. (|32l33l45p . The resonance is sufficiently nar- 
row {R* /\ai,f,\ ^ 34) for having a well defined interval of 
detuning where the intermediate regime of Eq. (j34p holds. 
This figure illustrates the fact that the effective range ap- 




-1,5 -1 -0,5 
(B-B^)/\AB\ 



FIG. 3: Spectrum of the shallow dimer for the narrow 
resonance of ^'^K at 752 G (abg = -35 ao, AB ^ -0.4 G, 
Sn — 1.5 fiB, Rvdw ~ 64.6 ao). Solid line: binding energy 
Edim obtained from the two-channel model in units of |5/iAZ3| 
as a function of the relative detuning {B — Bo)/\AB\; dashed 
line: approximate law for narrow resonances from Ea. (|32|l . 
Dotted line: approximation in Eq. (|33|l . derived in the in- 
termediate regime of narrow resonances; dashed dotted line: 
binding energy from Eq. (|45|l : dashed doubly dotted line: 
binding energy in the effective range approximation (|44p . 



proximation is relevant only for a small enough relative 
detuning \{B ~ Bo)/AB\ ^ 1) while the low energy ap- 
proximation in Eq. (j32p remains a good approximation 
of the exact solution of Ec^. (1^^ of the model even for a 
finite relative detuning. 



IV. THREE-BODY RECOMBINATION 

A. Wave equation for three bosons 

In this subsection, we derive the integral equation 
which permits to solve the general three-body problem 
of the present two-channel model. It is thus the first step 
in order to compute three-body recombination processes. 
We consider three bosons in their center of mass frame 
at energy E. The state of the system is a coherent super- 
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position of three atoms plus one atom and one molecule: 



I*) = 



d^Kd^k , t t t , 
d^K 



(27r)6 



/3KfoL«iKlO)- (47) 



(27r)3 



In Eq. (|T7l). /3k (respectively ^K.k) represents after 
symmetrization the atom-molecule wavefunction (respec- 
tively the wave function for the three atoms in the open 
channel, where K and k are two Jacobi coordinates in 
momentum space). In what follows, without loss of gen- 
erality ^K,k is considered to be an even function of k. 

Projection of the stationary Schrodinger's equation 
{E — -ff = 0, at energy E on the states with one atom 
and one molecule provides the equation for /3k ^ 

{E,,i - E^oi)Pk - 2A7K = (48) 

where for convenience, we introduced the function; 



7k 



d^k 
(2^ 



■Xk 



(^K,l 



2A 



hk,-^-J 



(49) 



In Eq. (|3H]), Ej-ci is the relative (or collisional) energy of 
a pair of center of mass momentum K: 



E„ 



E-- . 

4 m 



(50) 



It represents the energy of a pair in its center of mass 
frame while the third particle does not interact with 
an atom of this pair. The projection of the stationary 
Schrodinger's equation on the states with three atoms 
gives: 



(E'loi - 2ek)AK,k - Axk/3K - .goXk7K = 0. 



(51) 



In order to simplify subsequent equations, we introduce 
the effective wavefunction /3^: 



/3i 



eff 
K 



K 



2A2 



[Srei + Sfi{AB -B + Bo)] 



(52) 



where the factor can be also written in terms of the 



molecular energy by using Eq. ([2^ . From Eq. we 
obtain: 

iE,,i ~ 2ek)AK,k ~ Axk/3^^ = 0. (53) 

Similarly to the two-body case [compare with 
Eqs. (|14I15I) ]. the general solution of Eq. ([53)) can 
be written as: 



A 



K,k 



A 



(0) 
K,k 



eff 



Axk/3|? 



Sroi - 2ek + iO+ 



(54) 



where A^*^^ represents a possible incoming wave of three 
free atoms, and is thus an eigenstate of the kinetic energy 
in the center of mass frame. A^*^-* is non zero only if the 



energy E is non negative, in this case the second term 
in the right hand side of Eq. ([5^ represents an outgoing 
spherical scattered wave, thanks to the iO"*" term in the 
denominator. It is clear that for a negative energy, the 
factor between brackets in Eq. (|51l) cannot vanish and 
the prescription E ^ E + iO+ can be omitted. 

From Eqs. (|48l54p . we deduce a closed integral equa- 
tion verified by the effective atom-molecule wave func- 
tion: 

mexp{-m,E,cib'^/h'^) g 



A7rh^f{E,,i) 

d^k 



= 2 



/3k" + ^ 



+k 



(27r)3 ek + fiK + eK+k - - ^0^ 



where, 



(0) 

7k 



d^k 
(2^ 



Xk 



A 



(0) 
K,k 



2A 



(0) 

2 4 



(55) 



(56) 



Remarkably this equation has a form similar to the Skor- 
niakov Ter Martirosian (STM) equation ^Slj . In the limit 
where — > 0, the kernel in the integral of Eq. (1551) coin- 
cides with the one of the STM equation. For the more 
realistic situation where b is finite, the cutoff function Xk 
permits to avoid the Thomas collapse [l^l . A remarkable 
feature of the two-channel separable model used in this 
paper is that the two-body physics is entirely encapsu- 
lated in the diagonal part of the integral equation (|55|) 
only through the two-body scattering amplitude f{E). It 
has been shown recently that this property remains true 
in the four-body integral equation derived from the same 
model pj]. 



B. Three-body recombination constant 

We consider the recombination of the three bosons into 
a Feshbach dimer and one atom in the limit where there 
is no background dimer. Thus, the following inequalities 
hold: 



a > 



and 



< by/TT. 



(57) 



The derivation of the recombination rate and recombi- 
nation constant is obtained along the same lines as in 
Ref. [l^- We study the scattering process in the center 
of mass frame, we thus assume a vanishing total mo- 
mentum for the incoming wave. Moreover, we consider 
the limit of vanishing kinetic energy in the center of mass 
frame which is relevant for the Bose gas in the degenerate 
phase. We denote the wavevectors of the three incoming 
atoms by k°, kj and kg and we introduce the total and 
relative momentum Kg = k^ -I- kj and ko = (kj — k2)/2 
for the atomic pair 12. The incoming wave function can 
be written as 



1^(0)) 



dKdk 

"(2^ 



'^K,k"iK+k"iK-k" 



«Lk|0), (58) 
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where 



,(0) 



(2^)65(K-Ko) 



(5(k - ko) + 5(k + ko) 



(59) 



In the Umit of vanishingly smaU but positive energy we 
have: 



A 



(0) 
K,k 



(27r)6j(K)<5(k), 



(60) 



and we always assume in the following that E — Q. In 
this limit, the source term in Eq. ([55)1 which is given 
by Eq. ^ is 7^^ = 3(27r)3,5(K). The atom-molecule 
wavefunction /3k can be written as 



/3k = 3\/87ri?*(a - abg)(27r)35(K) 

A-Kh{K) 

+ VPcloscd^2_^2,^^_,0+- ^^^^ 

The coefficient in front of the delta distribution in 
Eq. ([5T|) is adjusted in order to cancel the source term of 
the generalized STM equation (|55|) . In the last term of 
Eq. (j6ip we introduced the wave number 



V3 



qdu 



(62) 



which corresponds to the relative momentum for the 
atom and the dimer formed in the inelastic three-body 
collision. Physically, the pole at i^T = i^dim + ^O"*" is 
linked to the outgoing spherical wave for the atom- 
dimer relative particle. The angular average appearing 
in Eq. ([55)) for the isotropic ansatz of Eq. (|61l) is per- 
formed in Appendix [A] After calculation one is left with 
the following inhomogeneous integral equation: 



/(i^rcl) 

00 

4 





"dim 



dk- 



Ki 



-C(i^,fc,0)e-i^'(^'+'=') 



247r?i^ ae~ 



VPcioscd mA 



K 



(63) 



where 



90HK) 
2A2 



3fi2^2 

4 m 



(64) 

We derive the recombination rate from the solution h{K) 
by considering formally the three-body system in a large 
fictitious box of size L. We consider the state of three 
particles in this box with periodic boundary conditions 
and defined by = |-0)/L^/^. The atom-molecule 

wavefunction in position space is obtained by a Fourier 
transform of Eq. : 



13 (r) = VPciosod47r 



(27r)3 



"dim 



(65) 



where r = r„ 



Tat and the dots indicate terms which 



are not relevant in the recombination process. This pro- 
vides: 



/3(r) = VPclosod h{Kdira)- 



(66) 



which corresponds as expected to an outgoing spherical 
wave for the relative particle composed of the atom and 
the molecule with the wave number i^dim- In the limit 
of large separation between the atom and the molecule, 
the outgoing wave function in the box is: 

^,\7(r)^V^^ ^^^';r^''7"" - (67) 

Hence, the probability current for the relative particle of 
reduced mass 2m/3 is: 



box^ 
out 



Aim 

J^closcd ^^dim 

L9 2m/3 



|/i(i^dim)lH- 



(68) 



The total flux associated with j is obtained by integration 
over the center of mass position and the full space angle 
(An) for the relative particle: 



L° m 



(69) 



The outgoing dimer has a probability Pcioscd to be in 
the closed channel, we thus divide this total fiux by the 
corresponding probability in order to get the probability 
flux of dimer formation. The atom loss rate is then 



jTvrbox 
"^'^dim 

dt 



= - — K^im\h{Kdim)\^^. (70) 



The general deflnition of the three-body recombination 
constant a-ccc for a system of N particles is: 



dt 



-"re 



N{N - 1){N 



(71) 



In the present case = 3 in Eq. (|7ip and by identiflca- 
tion with Eq. ([70]) we obtain the link between the recom- 
bination constant and the atom-molecule wave function: 



— Xdim|/l(i^dim)|' 

m 



(72) 



Finally, the determination of the recombination constant 
reduces to the calculation of h{Kdini); that is to the calcu- 
lation of the residue of the atom-molecule wavefunction 
at the wave number K = i^dim- 

We now apply Eq. (|72p to the case of narrow reso- 
nances in the regime R* S> |a| where the integral term 
in Eq. ([63]) is negligible for asymptotically large values 
of R* . We evaluate /i°*^(it'dim) to the lowest order by 
neglecting the integral term, thereby obtaining 



h'^\K) ^ ^U^^/(i?.el)(if^ - i^Jim 



)• (73) 
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Large values of the width radius R* are reached in the 
intermediary regime of Eq. ([M)) . so that we use the ap- 
proximate law Eq. (|33|) for the dimer's binding energy. 
Moreover, we use the approximation Pcioscd — 1 and we 
find the dominant contribution; 



dim J 



mA 



a(a - 



(74) 



Finally, in the intermediate regime the recombination 
constant can be approximated by the law: 



3R* (a — flbg)" 



(75) 



i?*3>|a| m 

This is the main result of this paper which is a general- 




0.01 



1 100 

R*/lal 



characteristics of the Efimov effect and was derived an- 
alytically in Ref. [1^ in the framework of a two-channel 
model strictly equivalent to the effective range approxi- 
mation. In Fig. ([5]), we plotted the three-body recombi- 




-10 
a (atomic units) 



FIG. 5: Three-body recombination constant a^cc in the region 
of small scattering length (|a| < |abg|) for the narrow reso- 
nance in ^^K in the vicinity of 752 G. Solid line: numerical 

(IZ3 



solution of Eq. (|63p : dashed line: approximate law of Eq. 
valid for large values of the ratio R* /\a\. 



nation constant for the same species in the region of zero 
crossing for the scattering length {B = Bq + AB) where 
it vanishes quadratically. 



FIG. 4: Three-body recombination constant arcc in units of 
h/m for the narrow resonance in '^^K in the vicinity of 752 G. 
Solid line: numerical solution of Eq. (|63[l for positive value of 
a; dashed line: numerical solution of Eg. (|63l) for abg < a < 0; 
dotted line: analytical solution of Ref. [2^ valid for large val- 
ues of the scattering length; dashed-dotted line: approximate 
law of Eq. (|75|) valid for large values of the ratio R* /\a\. In 
order to guide the eyes, this last curve is also plotted in the 
region where i?*/|a| is small. 

ization of Eq. (16) in Ref. 20] obtained in the effective 
range approximation and where the supplementary con- 
dition a ^ flbg was assumed. The deviation of Eq. ([75)) 
from this last law is substantial in the case where the 
background scattering length is negative and the ratio 
R*/\a\ can be arbitrarily small for a fixed value of R* 
while a shallow dimer still exists. In this regime, the re- 
combination constant tends to zero with an anomalous 
quadratic law in the scattering length a (the usual law 
brought by dimensional analysis is in a^). This regime 
can be reached for instance in the case of the narrow res- 
onance in potassium near 752 Gauss (see Fig. [3]). We 
illustrate this example in Fig. (U)), where the asymptotic 
law in Eq. (|75p is compared to the exact numerical cal- 
culation. In the region where the scattering length a 
is large (and positive) with respect to the width radius 
R* , we recover the oscillatory law in log-scale which is 



V. CONCLUSIONS 

In this paper, we used a separable two-channel model 
for the description of the magnetic Feshbach resonance in 
ultracold atoms. We studied the deviation of the dimer's 
binding energy with respect to the general threshold law 
Eq. ([28]) in the different regimes that can be achieved ex- 
perimentally. For narrow resonances, we defined a inter- 
mediate regime of detuning where the background scat- 
tering length is a relevant parameter. In this regime, we 
derived approximate laws for the binding energy of the 
shallow dimer and for the three-body recombination rate. 
These results show the interest of using a two-channel 
model for the description of ultracold gases in the reso- 
nant regime. In particular it appears as an useful tool for 
providing universal properties (that is, properties which 
are not explicitly dependent of the short range details 
of the interactions) in regimes where the effective range 
approximation is not quantitative. 
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Appendix A: Angular averages 

The ansatz for /3k used in this paper is isotropic. Thus 
in the homogeneous part of Eq. (j55p one has to perform 
the following angular average: 



47r q'^ + K^ + k'^ + K-k 



(Al) 



(in the recombination process considered in the present 
paper, q = 0). To perform the angular integration we use 
spherical polar coordinates of polar axis K/K. After the 
integration over the azimutal angle, we are left with 



I{K,k,q) 



C{K,k,q) 
Kk ' 



(A2) 



In Eq.dJU, C(K,k,q) is defined by: 
1 



j dx e~ 



— [Eiivt-t)-Ei{vt + t)] (A3) 



where Ei (z) = ds ^-j— is the exponential integral, 
and the variables {v,t) are: 



q^ + K^ + e 
Kk 



t = b^Kk. (A4) 



In the zero-range limit 6 — )■ we obtain the simple result: 



^ 1 , fq^ + K'^ + k'^ + Kk 
C = - m 



2 \q^ + K + k^ - Kk 



(A5) 



which coincides with the result of Skorniakov and Ter 
Martirosian ,31il. 
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